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Topological Anderson transitions, which are direct phase transitions between topologically distinct 
Anderson localised phases, allow for criticality in ID disordered systems. We analyse the statistical 
properties of an emsemble of critical wavefunctions at such transitions. We find that the local 
moments are strongly inhomogeneous, with significant amplification towards the edges of the system. 

In particular, we obtain an analytic expression for the spatial profile of the local moments which is 
valid at all topological Anderson transitions in ID, as we verify by direct comparison with numerical 
simulations of various lattice models. 
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The interplay of topology and disorder in electronic 
systems can give rise to critical transitions between lo¬ 
calised phases [IHS]- This is of particular interest in ID, 
as here there are good reasons to expect a system of elec¬ 
trons to realise an insulating phase, even in the absence 
of electronic correlations. On the one hand, a ID chain 
is generically unstable to dimerisation via Peierls distor¬ 
tion, while on the other, any amount of extensive dis¬ 
order leads to Anderson localisation. Remarkably how¬ 
ever, these two effects can be tuned in combination to in¬ 
duce critical behaviour. This occurs as distinct dimerisa¬ 
tion patterns realise distinct topological phases, whereas 
strong disorder washes out the detail of the dimerisa¬ 
tion, which drives the system between the phases. In 
general, such critical transitions occur in ID for disor¬ 
dered systems which realise symmetry-protected topolog¬ 
ical phases, and are described by two-parameter scaling 
[5]. We shall refer to the transitions as topological An¬ 
derson transitions, and it is the resulting critical states 
that are of interest to us in this paper. 

Whether a system may realise distinct topological 
phases is dictated by both its dimension and the anti¬ 
unitary symmetries of its Hamiltonian [7HS]. In partic¬ 
ular, for each dimension a certain five of the ten sym¬ 
metry classes of the Altland-Zirnbauer classification HO] 
allow for non-trivial topological phases. In ID these are 
the three chiral classes (AIII, BDI, CII) and two of the 
Bogoliubov-de Gennes classes (D, Dill). The situation in 
ID is however more special as here, irrespective of sym¬ 
metry class, all topological Anderson transitions share 
the same critical properties, a phenomenon dubbed supe¬ 
runiversality |3]. This was established by identifying ex¬ 
act mappings between representative single-channel mod¬ 
els for each of the five classes, and by showing that for a 
long chain the properties of a system with many channels 
are determined by a dominant effective single channel. 

It is of great interest to understand the properties 
of the critical wavefuntion at this superuniversal tran¬ 
sition. This has been partially investigated by Balents 
and Fisher m who calculated the system-size scaling of 
the moments of the normalised wavefunction in the bulk 


and their spatial correlations 
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where (...) denotes an averaging over a disorder ensem¬ 
ble, N is the system size and g > 0. The 1/Wscaling 
of the moments is characteristic of localisation, while the 
algebraic decay of the spatial correlations reflects the crit¬ 
ical nature of the wavefunction |I2j . 

Here we analytically show that the local moments are 
in fact inhomogeneous throughout the system. In the 
bulk they scale as 
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and directly at the edges they scale as I/'/N. Equation 
(U describes how the 1/A^-scaling in the bulk changes 
to I/VN towards the edges. The strong amplification 
of the disorder-averaged wavefunction near the bound¬ 
aries (|/(l)/^)/(|i/;(Af/2)/'J) oc v/V distinguishes it sig¬ 
nificantly from all other critical wavefunctions, and can 
be understood as a remnant of the topological edge mode 
that is created/destroyed by the transition. 

These results are derived from a specific single-channel 
model in class BDI for which the critical wavefunction at 
the transition is exact and analytically tractable. Based 
on superuniversality, we argue that the results are gener¬ 
ally valid at all ID topological Anderson transitions. We 
verify this by directly comparing them with numerics for 
both a class BDI model with no A = 0 wavefunctions 
away from the N ^ oo limit, and for a multi-channel 
model belonging to symmetry class D. 

To explore the critical physics of ID topological Ander¬ 
son transitions we focus on a disordered dimerized tight- 
binding chain. This model appears as the static limit of 
the Su-Schrieffer-Heeger model mui], which is used for 
example to describe polymers such as polyacetylene. For 
reasons that will be made clear shortly, we primarily fo¬ 
cus on a chain of odd length 2N -\- 1, and label the sites 
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from 0 to 2N\ 


the normalized zero-energy wavefunction: 


N 

H = [ — tjc\j_ 2 C 2 j-l — tjc\j_iC 2 j + h.c.] . (3) 

i=i 

Disorder appears through the hopping amplitudes 

tj = t + Wujj , tj =t' + W'ujj , (4) 

where ujj , w' are independent random variables uniformly 
distributed on the interval [—0.5, 0.5]. The model is in¬ 
variant under the time reversal symmetry T = K, where 
K is complex conjugation, and has a chiral symmetry 
SHS~^ = —H, with S = parity op¬ 

erator P = TS squares to the identity when acting on 
single-particle states, and so the model belongs to the 
BDI symmetry class. 

The odd length of the chain guarantees the existence 
of a zero-energy eigenstate for any set of parameters tj, 
t', [la [16]. This follows from the chiral symmetry of 
the model which restricts the spectrum to be symmetric 
about zero. The Schrodinger equation at zero energy 
reduces to a simple recursion relation with solution 

i’2j = V'O , '02j-i = 0, j = 1,..., iV, (5) 

where the ^ indicates that this wavefunction is unnor¬ 
malised. This eigenstate describes a topological bound¬ 
ary mode that is localised at one edge of the chain. For 
the clean case (tfe = t, tJ,, = t') the mode resides at the 
left edge if jtj < jt'j and at the right edge if jtj > jt'j, while 
precisely at jtj = jt'j there is a topological transition at 
which the wavefunction is uniform. In the presence of 
disorder this transition becomes a topological Anderson 
transition, and it can be driven purely by a change in the 
disorder strength. 

It is worthwhile to contrast this behaviour of the edge 
mode with that for an even length chain, for which there 
are boundary modes at both edges of the chain on one 
side of the transition, and no boundary modes on the 
other side. The connection can be made by considering 
adding an extra site to the odd length chain. On its own 
the extra site has a single localised zero energy mode. On 
coupling it to an end of the chain this localised mode will 
either remain (protected by the chiral symmetry) if there 
is no zero energy state at that edge, or it will combine 
with and annihilate the boundary mode if there is one 
present. More formally, for the even length chain the 
two sides of the transition can be distinguished by the 
Zak-Berry phase [imiH], which is a Z 2 bulk invariant, 
and the extra site (which breaks the 2 -site unit cell) can 
be thought of as giving rise to a defect mode [19]. We 
remark that for the odd length chain there are two choices 
for the unit cell, and so there is no natural way to call 
one phase trivial and the other topological. 

We now calculate the local moments per unit cell of 


leAn 

where n indexes the unit cells of the underlying lat¬ 
tice, and A„ is the set of sites in the n**' unit cell. 
In the present case Iq{n) = (|' 02 nP®), with ip2n = 

' 02 ra/(X)m=o l^ 2 mP)^^^, as the wavefunction has support 
on only the even sites of the chain. The normalisation 
factor contains non-trivial dependence on the disorder 
parameters , ui'f., and to handle this we use a technique 
developed for continuum models that maps the calcula¬ 
tion to one in Liouville quantum mechanics [IIIEO]. In 
the present case, working explicitly with a lattice model 
allows for a controlled description of the problem in the 
strong disorder regime. We outline here only the main 
steps of the calculation and defer the details to the Ap¬ 
pendix. 

The disorder averaged local moments of the wave- 
function Iq{n) can be expressed as a functional inte¬ 
gral over random walk configurations of the variables 
Zk = 2 In In the thermodynamic limit TV —> cx), 

the behavior of this integral is determined by the first 
two moments of the distribution function of the Zk- As a 
result, evaluating Iq(ji) reduces to solving a discrete-time 
version of Liouville quantum mechanics - the quantum 
mechanics of a particle moving in an exponential poten¬ 
tial. This latter problem can be solved by approximating 
the exponential potential by a hard wall. 

As shown in the Appendix, the averaged local moments 
can be expressed in terms of a propagator P{x, x', m) as 


Iq{n)=AqJ J J dxdydzP{z,y,n)V{y)P{y,x,N - n), 
000 

( 7 ) 

where Aq is a normalisation constant and V(x) = 
g- 2 gV^x ^ ^ (^^ 2 ^ _ Within the hard 

wall approximation, the propagator P{x,x',m) is given 
by a solution of the classical drift-diffusion equation with 
an absorbing boundary condition at cc = 0 : 

P{x,x ,m) =— e ™ )e ™ , ( 8 ) 

V7rm 


which satisfies the initial condition P{x, x', 0) = 5{x — x') 
for X, x' > 0. A non-zero value of the drift constant 
7 = — (zk) I2yja corresponds to the insulating phases, 
which are distinguished by the sign of 7 , while zero drift 
corresponds to criticality, and thus {zk) = 0 is the con¬ 
dition on the distribution of tj, t'j that the wavefunction 
of Eq. ([^ is critical. 

We briefly summarise the behaviour of the moments 
away from the critical regime. Evaluating the integrals 
in Eq. Q for 7 ^ 0, we obtain 


Iq{n) 


C' ' p-^V^qjn ^ > 0 
qTiq)^ > 
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where C = 1 — e , F is the gamma function, and we 
assume that 'yn ^ 1 for 7 > 0 , and | 7 |(iV — n) 

^/N — n ^ 1 for 7 < 0. Therefore the wavefunction is ex¬ 
ponentially localized at the left edge for 7 > 0 and at the 
right edge for 7 < 0. Thus the two topologically distinct 
insulating phases in the model can be distinguished by 
probing the local statistics of the wavefunctions. In terms 
of the underlying diffusion process, the spatial asymme¬ 
try of the wavefunctions can be interpreted through the 
interplay of the non-zero drift with the absorbing bound¬ 
ary condition at a; = 0, see the Appendix. 

For critical wavefunctions, corresponding to 7 = 0, 
integrating Eq. 0 yields 


Iq{n) 


Tr2i ^q^T{q) 


n, N — n 1. (10) 


The scaling of the moments is Iq{n) (x. 1/N deep in the 
bulk of the system, where n oc. N, which is the scaling 
law typical for a localized wavefunction |12j . At the same 
time, the moments get strongly amplified approaching 
the edges of the system. Exactly at the edge, the mo¬ 
ments are given by 


^ 9 ^^) 2<!q2r{q)y/^y/N' 

The 1/A^'^-scaling of lq(0) with a non-integer exponent 
T = 1/2 is the scaling law specific for fractal wavefunc¬ 
tions [12]. The fact that r < 1 is however unusual, as 
typically amplitudes of the critical wavefunctions are sup¬ 
pressed rather than enhanced near the boundaries and 
T > 1 jnj. The inhomogeneous nature of the local mo¬ 
ments and the amplification of the critical wavefunction 
at the edges can be traced to the presence of the local¬ 
ized boundary modes in the insulating phase. In terms of 
the classical diffusion process it can be traced to the role 
of the potential V{x), which gives an increased weight 
to trajectories passing near the absorbing boundary con¬ 
dition at a; = 0. Viewed in this way it is natural that 
Iq{0)‘^/Iq{N/2) oc see the Appendix. 

The spatial correlations of the critical wavefunctions 
can be characterised by the two-point correlation func¬ 
tion Cmn = (|'02m'02n|®)- Using the same approach as 
above, we obtain that Cmn oc |to — in the 

bulk of the system, and Con oc at the edge. 

The power-law scaling of the correlation function is an¬ 
other manifestation of its fractal nature HD. 

To investigate the range of validity of the above results, 
which are derived in the limit of large V, we analyse the 
critical wavefunction numerically. For this we examine 
a critical point of the Hamiltonian of Eq. ([^ given by 
t = 1, t' = 0, W = 2, W' = 4, and compute the local 
moments by averaging the zero-energy wavefunction over 
10® —10® disorder realisations. The system-size scaling of 
the local moments, both directly at the edge and in the 
centre of the chain, are presented in Fig.j^ It is seen that 
the edge scaling is more robust than the bulk scaling for 



FIG. 1: (Color online) The system-size scaling of the local 
moments Ig(n) of the critical wavefunction of Eq. © com¬ 
puted (i) at the edge and (ii) at the centre of the chain. The 
q = 1/4,1/2,1, 2,4 moments are plotted, and the dashed lines 
give fits to (i) l/vfv scaling at the edge and (ii) 1/A scaling 
in the bulk. 
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FIG. 2: (Color online) The spatial profile (n dependence) of 
the local moments of the critical wavefunction of Eq. © with 
q = 1/4,1/2,1, 2,4 plotted for (i) A = 50, (ii) A = 100, (hi) 
A = 200. The dashed lines show fits of the data to Eq. ( |10[ ) 
where the normalisation is set by fixing the value at the centre 
of the chain. 


small system sizes, while for increasing JV an excellent fit 
to Eq. ( jTTj ) and ( 101 , for the edge and bulk respectively, 
emerges. 

In Fig. the spatial prohle of the q = 1/4,1/2,1, 2,4 
local moments for A = 50,100,200 are plotted, and com¬ 
pared with a fit to Eq. (10). It is seen that for g > 1 
the formula captures the moments right out to the edge. 
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FIG. 3: (Color online) Plots of the local moments for the 
even-length chain, obtained by setting = 0 in Eq. (|^, with 
the disorder average taken over wavefunctions whose energy 
lies within a window about 0, as described in the text. Plot 
(i) shows the system-size scaling of the edge moments for q = 
1/4,1/2,1,2,4 and the dashed lines give fits to 1/v/V scaling. 
The data is terminated at the black squares on the g = 1, 2,4 
lines as beyond these system sizes numerical inaccuracy arises 
due to machine precision. In (ii) the spatial profile of the 
q ~ 1 local moments are plotted for N = 25, 50,100, 200 and 
the dashed lines show fits of the data to Eq. ( |10[ ) where the 
normalisation is set by fixing the value at the centre of the 
chain. 


while for g < 1 there are deviations near the edges of the 
chain which get smaller as the system size is increased. 
Since the moments become uniform in the limit g —>■ 0 
one would need diverging system sizes in order to see the 
predicted form in this regime. 

We now perform checks on the generality of our an¬ 
alytic expressions, to support the claim that they are 
valid at all topological Anderson transition in ID. First 
we consider a dimerised tight-binding chain with an even 
number of sites, as then there is no wavefunction exactly 
at if = 0 for any hnite length. To examine the criti¬ 
cal properties in this case we take the disorder average 
(...) to run over all wavefunctions in the disorder en¬ 
semble whose energy is within a window [—Fltyp, i?typ], 
where Atyp = exp((ln(minQ |ifa|))) with {Ea} the set 
of eigenenergies for a given realisation. The even length 
chain is realised by setting ti = 0 in Eq. ([^, which decou¬ 
ples the site 0, and it is analysed using exact diagonalisa- 
tion. In Fig. I^i) the system-size scaling of the local mo¬ 
ments at the edge of the chain are plotted and compared 
with the predicted 1/\/N scaling, while in Fig. I^ii) the 
spatial profile of the <7 = 1 local moments are presented 
and compared with the prediction of Eq. ( [Io| . In both 
cases a very good fit is seen, which indicates that the re¬ 
sults are valid in general and not just to exactly critical 
wavefunctions. 


FIG. 4: (Golor online) Plots of the local moments for the 
Hamiltonian Hd with to = 1/4 and Ach = 5. The quantities 
plotted are the same as those in Fig.[^ and identical behaviour 
is observed. 


The second check on the generality of our results that 
we perform is to examine a multi-channel model in the 
Bogoliubov-de Gennes symmetry class D. Specifically we 
consider the following Hamiltonian 

2N 

- i c]B,C, + h.c.] , (12) 

J=o 


where fu is a hopping parameter, Cj is a column vec¬ 
tor {cj,i, Cj_ 2 ) • ■ •) Cj^AT,,^} of Nch fermions at site j, IE 
is an Nch x Ach random matrix satisfying Bj = —Bj, 
and Nch is the number of channels. In practice we set 
Bj = Aj — Aj where Aj are real matrices whose ele¬ 
ments are independent random variables uniformly dis¬ 
tributed on the interval [—0.5,0.5]. The Hamiltonian be¬ 
longs to symmetry class D when A/h > 2, as complex 
conjugation K provides an anti-unitary symmetry that 
anti-commutes with the Hamiltonian, whereas in the ab¬ 
sence of a sublattice symmetry there is generically no 
anti-unitary symmetry that commutes with the Hamilto¬ 
nian. The model sits exactly at a topological Anderson 
transition for any tjj, and can be taken into a topological 
phase by adding a dimerisation to the hopping. Indeed, 
to view the topological properties of the model one must 
regard it as having a two-site unit cell, and it is for this 
reason that we set the sum to run to 2N in Eq. (12|. 

Again we numerically examine the disorder averaged 
local moments of the critical wavefunction. For this 
we take Nch = 5 and obtain the critical wavefunctions 
through exact diagonalisation. In Fig. we plot again 
the system-size scaling of the local moments at the edge 
of the chain and the spatial profile of the < 7=1 local mo¬ 
ments, as were presented in Fig. We find once more 
an excellent ht, and view this as confirmation that our 
results are generally applicable at topological Anderson 
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transition in ID. 

In summary, from the known exact zero-energy wave- 
function of the Su-Schrieffer-Heeger model we have de¬ 
rived analytical expressions for the system-size scaling of 
the local moments, which are valid at all topological An¬ 
derson transitions in ID. The generality of the expres¬ 
sions follow from the superuniversality of these transi¬ 
tions, and is also evidenced by our numerical simulations 
of other models, which show excellent agreement despite 
having distinct features from the original model. The 
results show that the moments are inhomogneous, with 
strong amplification towards the edges of the system, and 


this distinguishes critical wavefunctions at topological 
Anderson transitions in ID from all other critical wave- 
functions studied so far. It is an interesting open problem 
whether similar features of critical wavefunctions can be 
found in higher dimensions. 
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Appendix A 

In this Appendix we provide details of the analytic 
calculations leading to the results quoted in the main 
text. 


1. Moments of the normalized wavefunctions 

Moments of the normalized wavefunctions and classical 
diffusion 


The explicit form of the wavefunction given by Eq. 
(5) allows us to straightforwardly calculate its statisti¬ 
cal properties, if the normalisation condition is not taken 
into account. However, for the normalized wavefunction 
the problem becomes much more complicated, as the nor¬ 
malisation constant depends on all variables tk, t'f. in a 
non-trivial way. In order to overcome this problem we 
map the calculation to one in Liouville quantum mechan¬ 
ics |AI1 IA2) . 

First, we define the random variables Xn = 
nfe=i we set a;o = 1 . Then the disor¬ 

der averaged local moments of the normalized wave- 
function Iq{n) = (|V’2nP'^) are given by Iq{n) = 

(^^n/iJ2p=o The denominator of Iq{n) can be rep¬ 

resented as 

/N \-9 oo 

, (Al) 

where T is the gamma function. With the help of the 
above identity we obtain 


1 

Iq{n) = —— / daa'^~^e~°‘J{a), 

^ ( 9 ) Jo 


(A2) 


where J{a) = (xfe “ ^p=i ^ and we used the fact that 
xq = 1. The disorder averaging can be performed using 
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the statistical independence of the random variable Zk = 


J{a)= / dzi ... dzN Pz{zi)... “^p= 


(A3) 


where Pz{zk) is the distribution function of Zk- The structure of the above integrand suggests to introduce new 
variables 


p 

^ Zk, Zp = tp- tp-i, for = 2,..., and zi = ti, (A4) 

k^l 

in terms of which J{a) takes the form 


J{a) 



N 

■■dtN WPz{tp 

P^l 




(AS) 


where to = 0, gp{t, a) = e for n and gn{t, a) = Using the expression for Pz{z) in terms of its Fourier 

transform F{s) = dze^^^P^iz), we obtain 


J{a) 



dsN 



N 

dti... dtN n a), 

p^l 


One can notice that 


^ 2Sp(ip tp — 2^) 


Fisp) 


< |F(sp)| <1 = F(0), 


(A 6 ) 


(A7) 


hence the main contribution to the integral in the limit N ^ oo comes from Sp —)■ 0. Therefore we can use the Taylor 
expansion of F{s) about s = 0. Keeping only the terms up to the quadratic order we can write F{s) « ^-^bs-as^ ^ 
where b = — (zk), a = {{zi) — {zk)‘^)(2. Then the integrals over Sp become Gaussian and can be evaluated explicitly: 


J{a) 




(tp —tp I +b) 


-\-OLe''P 


(AS) 


This expression can be viewed as a discrete in time path integral for a particle moving in an exponential potential - 
Liouville quantum mechanics. One can approximate such a potential by a hard wall potential: 

e-“^‘ « 9{z{a) - t), (A9) 

where 9{t) is the Heaviside step function, = C and C is a constant of order 1. Within this approximation J{a) 

is given by 


J{a) 



g9ir.g-Ep 


(tp—tp I +b)‘^ 



Vtt 


(AlO) 


where Vp = {z — tp)j2y/a^ Pq = zl2^~a and 7 = bj2y/a. In order to calculate the above integral, we introduce an 
integral operator K defined as 


{Ku){x) = 9{x) 




u{y)- 


The propagator corresponding to the path integral (AlO) is related to K as 


P{x,x',m) = (x), 5x'{x) = S{x - x'). 


(All) 


(A12) 
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If P{x,x',m) is known, then J(a) can be found as 

pOO poo 

= dxi dx2P{ro,X2,l2)V{x2)Pix2,xi,li), 

Jo Jo 


(A13) 


where V{x) = e ^ _ iy_— n, I 2 = n. In the hard wall approximation P{x,x',m) = 0 for x < 0, and hence 

the integration over a in Eq. (A2) must be restricted to a < C. Taking into account that = C^ja, the 


expression for the moments reads 

Cq 

/,(") = 


da - / dxi 

Oi. ./n 


dx2P{ro, X2, kW (X2)P(X2, Xi, Zi). 


(AM) 


Changing the integration variable a by z{a) and noting that e « 1 in the hard wall approximation for 2 ; > 0, 

we obtain 


Iqin) 


2^0(79 

m 




dxsP^Xs, X2, ^2)V^ {X2)P{X2,XlJl 


(A15) 


where we finally changed z by X 3 = xq = zj^^fa. This equation is equivalent to Eq. (7) in the main text with 
A, = 2v^C'Vr(Q). 

In order to find the propagator P(x,x',m), we notice that, as follows from its definition, it satisfies an integral 
equation 

P(x,x',m + I) = d(x) [ P{y^x',m). (A-16) 

Jo 'J'x 

Assuming that P{y,x',m) is a slowly varying function compared to (this assumption is justified by the 

solution that we obtain, so is self-consistent), one can transform the above integral equation into the drift-diffusion 
equation with the absorbing boundary condition at x = 0. The solution of the drift-diffusion equation obtained by 
the standard eigenfunction method reads 


P(x, x', m) = 


{x — x' + 'ym)‘^ 


{x — x'-\-~^m')‘^-\-Ax 


— e 


(A17) 


One can check by direct calculation that the above expression for P{x,x',m) solves the integral equation (A16) 
provided that m,x^ 1, i.e. it works well for large times and far away from the boundary x = 0. Eq. (AI5| and Eq. 


(A17) give a complete description of the local moments of the wavefunctions in terms of classical diffusion. 


b. Localized wavefunctions 


The insulating phase is characterized by a non-zero drift constant 7 7 ^ 0. Substituting Eq. (A17) into Eq. (A15) 
and integrating over xi and X 3 we obtain 

Iq{n) = f dx (1 -I- erf (u+(x)) — e“‘^'^'’'erfc (u_(x))) (1 -f erf (ti_(x)) — e‘^'^'’'erfc (t>+(x))) (A18) 

2f (Q) Jo 

where u±(x) = (x ± 7li)/Vk, x±(x) = (x ± 'yk)/^^, erf (x) is the error function and erfc(x) = 1 — erf (x). To 
evaluate this integral we may use the following approximation for erf (x) 


erf (x) « 26{x) — I, (A.19) 

which is justified in the above integral provided that "fli ^ ^ I. This approximation yields the following result 

for Iq{n) 


Iq{ri) 


CJ -2y/dq-in a/ > 0 

qT(q)^ ’ 

C" g-2%/a<;|7|(Af-ri) < Q 

qT{q) ’ ' 


(A20) 


The positive 7 case corresponds to the wave function localized at the left edge, while the negative 7 case corresponds 
to the wave function localized at the right edge. Thus the two different topological phases in our model can be 
distinguished by probing the local statistics of the wavefunctions. 
















In order to determine C we use the normalization condition 




1 , which gives 


C = 1 - 


(A21) 


The spatial asymmetry of the result can be naturally understood in terms of the underlying diffusion process. Let us 
consider for simplicity 7 > 0 and two sub-cases n = 0 and n = N. For positive 7 , the propagator P{x, x', m) de scribes 
trajectories with a drift to the left (the mean value of x decreases in ’’time” m), as it follows from Eq. (A17). Now 
consider the wavefunction at n = 0 , the moments are given by 




2 v^C'« 

m 




dX2V{X2)P{X2, Xi,N). 


(A22) 


The propagator P{x 2 ,xi, N) describes trajectories of classical particles, which start from x = xi a.t time m = 0 
and propagate to a; = a ;2 at m = A^. The contribution of each trajectory is weighted with an exponential function 
^( 3 ^ 2 ) = at the final point. For this reason the main contribution to the results comes from the trajectories 

with X 2 close to zero. Such trajectories can be found and, due to the drift term, the most probable starting point for 
them is xi = 'yN. The contribution of such trajectories yields the result of order of one. 

If we now consider the opposite edge n = N, then 

4(-^) = -p/ N / dxi / dx 2 P{x 2 ,Xi,N)V{xi). (A23) 

r(g) Jo Jo 


The weight function V{xi) = jg associated with the starting point, which forces the initial coordinate xi to 

be close to zero. On the other hand, the absorbing boundary condition at a; = 0, together with the presence of the 
drift term which tends to decrease the position of a particle in time, lead to the exponential small survival probability 
of such trajectories and to the exponentially small final result. 


c. Critical wavefunctions 


The critical phase is characterised by a vanishing drift constant 7 = 0 . In this case integration over x^ and Xi in 
Eq. (A15l leads to the following result: 


Iq{n) = 


2 v^C'« 

m 


dx erf 




6 rf 

Vve; Vve 


-2qy/ax 


(A24) 


Recalling that li = N — n, I 2 = n and assuming that the site n is in the bulk of the system, we are allowed to replace 
two error functions in the thermodynamic limit iV —>■ 00 by their asymptotic expressions for li,l 2 —>■ 00 


erf( 4 ) = ^+0 

.VlJ VtW 


43 / 2 ) . 


Then the remaining integration over x yields 

Ig{n) = 


2 C« 


q^T[q)TTayJn{N — n) 


(A25) 


(A26) 


In order to fix the value of the constant C we use again the normalization condition = I, which gives 

C = al2. 

If n = /3N with A^-independent constant /3, then Iq{n) oc 1/N, which is the scaling law typical for a localized 
wavefunction. 

However, if we are interested in the scaling of the local moments on a site which is close to the boundary, the result 
will be very different. Indeed, let us consider n = 0. Then using Eq. (A22) and the same steps as above, we obtain 


4(0) = 


Cq 


q^T{q)ynTd\/N 


(A27) 


Thus 4(0) oc \/y/N, which is the scaling law typical for a fractal wave function. The same scaling is valid for the 
opposite edge n = N. 
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The two different results for the scaling of the moments at the edge and in the bulk of the system can be traced 
to the presence of the localized boundary modes in the insulating phase. It can be also understood in terms of the 
classical diffusion process as follows. The propagator P{x2,xi, N) describes traje ctories of classical particles, which 
start from x = xi at time m = 0 and propagate to x = X2 at m = N. In Eq. (A22 1 the contribution of each trajectory 
is weighted with an exponential function V{x 2 ) = at the final point . For this reason the main contribution 

to the result comes from the trajectories with X2 close to zero. Thus Eq. (A22) gives a probability of getting from any 
starting point Xi to an X 2 in the vicinity of zero in time N, and according to our calculations this scales as 
Similarly, Eq. (A15l for n = N/2 counts the contributions from the trajectories starting at any point xi reaching an 
X 2 in the vicinity of zero in time N/2, and then reaching an arbitrary final point X 3 in time N/2. The probability for 
the first part of the process scales as 1/\/N, but the same is true for the second part of the process due to the time 
reversal invariance of the problem. Hence the total probability must scale as {1/'/N)^ = 1/N. 


2. Two-point correlation function 


Spatial correlations of the wavefunctions can be determined through the two-point correlation function 

Cmn = • 


(A28) 


Performing all the steps which led us to Eq. (A15) we obtain 
2 ^C‘i 


c — 

'-^m.ri. — 


m 


pOO poo pOO pOO 

/ dxi / dx 2 / dx 3 / dXiP{Xi,X 3 , 13)11 {x 3 )P{x 3 ,X 2 ,l 2 )U{x 2 )P{x 2 ,Xi,li), (A29) 

Jo Jo Jo Jo 


where U{x) = e ^ li = N — n, I 2 = n — m, I 3 = m. Integration over xi and X 4 , and expansion of the resulting 
error functions gives 


C'mn — 


T^^{q)Vhh Jo 


poo poo 

/ dx2 / dX3X3e~‘^'^'^^P{X3,X2,l2)x2e~'^'^'^'^. 

Jo Jo 


(A30) 


The above integral can be computed explicitly in terms of the error functions. However, in order to find how it scales 
with |m — n|, it is more convenient to write it in the momentum space using the complete set of standing waves 


{z\k) — y ~ sm(A:z), fc > 0 . 

Writing the propagator in the momentum representation 

pOO pOO 

P(x3,X2,l2) = {x3\Pih)\x2) = dk dfc'| fc) (/c |P(l2) | fc') (fc'^ 2 ), 

Jo Jo 

calculating the matrix elements 

poo poo 

{k\P{l 2 )\k') = — / dx dys\n.{kx)s\n{k'y)P{x,y,l 2 ) = ^'^5(k — k'), 

Jo Jo 


and substituting two equations above into Eq. (A30), we find 

SVaC? 


Cm.n, — 


q poo / POO 

= / dke-^^’^ P / dx{x\k)xe-<i'P^^ 

1^3 Jo \Jo 


'n-^iq)Vhk . 

Evaluating the integral over x and scaling the variable k —>■ qy/ak we obtain the final result 


n — 


i-oo j^ 2 g-i 2 <? 

Tr'^aq^T{q)y/I/T/} Jq + I)"* 


(A31) 


(A32) 


(A33) 


(A34) 


(A35) 


which coincides up to a constant with the result by Balents and Fisher |A2) . For I 2 = 0, the integral over k yields 
7r/32 |A3] and hence we reproduce the expression for the moments (A26). For I 2 ^ 1 and m and n of order of IV ^ 1 
we find 


Cmn oc |m - n| \ 


(A36) 
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showing that the critical wavefunctions exhibit power law correlations in the bulk of the system. 

If o ne of the points is located at the edge of the system m = 0, then the above derivation should be modified. Eq. 
(|A29|) must be replaced by 


Con — 


2 y/aC'i 


dxi / dX2 / dX3U{x3)P{x3,X2,l2)U{x2)P{x2,Xi,li), 


where li = N — n, I 2 = n. Performing all the steps similar to the ones, which led to Eq. (A34), we find 


Cqu — 


poo poo poo 

= / dke~^^^ / dx3{x3\li)e~’^'^^^ / dx2(x2\k)x2er’^'^^'^. 

1 do do do 


Jo 

Computing the integrals over X 2 and X 3 we arrive at the final result 

lec? 


Con — 


f°° ofc /4 

Tr^/^y/dq‘^T{q)^/]J do ^ + 1)^ 


(A37) 


(A38) 


(A39) 


For I 2 = 0, the integral over k gives tt/IG and we reproduce Eq. (A27). For Z 2 S> 1 and N — n oi order of N we find 


Con OC 

which has the same scaling with m as in the bulk, but a different scaling with N. 


(A40) 


Finally, for the both points located at the opposite edges an analogue of Eq. (A29) reads 


2 faC^ 

Con = -p. ^ / dx2 / dx3U{x3)P{x3,X2,N)U{x2), 

^ \Q) Jo Jo 


which can be transformed to the momentum representation as 


Con = 


2v^C« 


dke-^^"/^ 


dx{x\k)e 


— qy/ax 


Calculating the integral over a;, we obtain 


Con — 


4C"^ 1 ^ 2 ^-Nq‘^aC /4 


rqTiq) 


dk 


(fc2 + 1)2 


(A41) 


(A42) 


(A43) 


which implies that Con N 
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